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SumrBjry 

The  timB-harmonic  excitation  is  either  an  interior  point  source  of 
sound  or  a  mechanical  point  force  located  on  the  pipe’s  surface.  The 
acoustic  power  radiated  into  the  exterior  fluid,  the  vibration  response  of 
the  pipe’s  wall,  and  the  pressures  in  the  exterior  and  interior  fluids  are 
computed  by  a  simple  integration  scheme.  Numerical  results  are  presented 
for  the  case  of  a  water'-filled  steel  pipe  that  is  surrounded  by  air. 
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■  cartesian  coordinates 

=  cylindrical  coordinates 
=  spherical  coordinates 

=  axial,  tangential  and  radial  displacements  at 
shell  mid-surface 

■  axial  and  circumferential  harmonic  wavenumbers 
=  transforms  of  shell  displacements 

=  interior  and  exterior  pressures 

=  radian  frequency  of  vibration 

=  density  and  sound  velocity  of  interior  fluid 

=  density  and  sound  velocity  of  exterior  fluid 

=  wavenumber  in  fluid,  w/c. 

1 

=  radius,  thickness  and  density  of  pipe 

■  Eh/Cl-o2).  E  is  Young’s  modulus  and  a  is 

Poisson's  ratio 

*  h2/12az 

=  Dirac  delta  function 

=  1  when  n=0,  and  2  otherwise 

=  amplitude  of  point  forcB  located  at 
cylindrical  coordinates  (a, 0,0) 

=  amplitude  of  point  source  located  at 
cylindrical  coordinates  (Xo,0,0) 

=  (X-Xq)2+Y2+Z2  -  r2*  X2  -  2rXg  cos(0)  +  Z2 

=  Bessel  functions  and  their  derivatives.  H  =J  +iY 

n  n  n 

=  (k|  -  a2)^  with  Im. (y^)  >  0  to  satisfy  the 
radiation  condition 

*  hysteretic  loss  factor  of  shell 
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INTRODUCTION 


A  previous  report  [  1  ]  has  analysed  the  theoretical  problem  of  far-f ield 
sound  radiation  from  fluid-filled  pipes  that  are  excited  by  point  forces  and 
point  acoustic  sources.  The  salient  points  that  emerge  from  the  work  are 
that  the  frequencies  at  which  'peaks’  appear  in  the  sound  spectra  depend  on 
the  location  of  the  observation  point,  and  that  the  far-field  radiation 
characteristics  are  understood  better  when  wavenumber  versus  frequency 
dispersion  plots  are  available.  In  practice,  however,  measurements  in  the 
far-field  are  not  usually  possible  because  the  environment  is  inevitably 
reverberant.  It  is  thus  necessary  to  average  spatially  the  pressure 
measurements,  in  a  suitable  way,  to  give  an  estimate  of  the  radiated 
acoustic  power.  Other  responses  of  central  interest  to  an  experimentally 
based  programme  are  the  pressure  fields  near  to  the  pipe's  surface  and  the 
acceleration  response  of  the  pipe’s  wall. 

In  Section  2,  the  pressure  fields  and  the  pipe's  displacements  are 
represented  by  Fourier  transforms,  and  the  acoustic  power  in  tne  surrounding 
fluid  is  obtained  from  the  stationary  phase  approximation  to  the  far-field 
pressure.  The  integral  representations  of  the  fluid  pressures,  pipe’s  dis¬ 
placements  and  radiated  power  are  evaluated  by  a  simple  numerical  quadrature 
scheme.  In  Section  3,  numerical  results  are  discussed  for  the  case  of  a 
water-filled  steel  pipe  that  is  radiating  into  the  surrounding  air.  The 
numerical  results  are  in  the  form  of  decibel  level  versus  frequency  plots 
for  the  responses  in  the  individual  circumferential  harmonics,  n=0,  1  and  2. 

The  numerical  results  are  better  understood  by  inspection  of  wavenumber 
versus  frequency  plots,  but  further  numerical  work  is  necessary.  First,  it 
is  necessary  to  search  the  literature  for  measurements  that  have  been 
obtained  under  controlled  conditions,  because  a  comparison  with  theoretical 
estimates  would  be  particularly  valuable.  Secondly,  plots  of  acoustic 
intensity  vectors  [  6  ]  would  help  to  illustrate  the  essential  physics  of  the 
fluid-pipe  system.  Thirdly,  the  eigenvectors  should  help  to  explain  the 
curious  maxima  and  minima  that  appear  in  the  spectra  of  the  responses  at  a 
distance  from  the  source,  Forthly,  the  range  of  applicability  of  the 
concept  of  matching  wavenumbers  of  in-vacuo  pipe  and  rigid-walled  duct  needs 
to  be  investigated  in  some  detail  because  it  is  being  used  increasingly  to 
explain  certain  aspects  of  measured  data.  Finally,  the  effect  of  fluid  flow 
on  the  stability,  vibration  and  sound  radiation  needs  to  be  investigated. 

2.  PROBLEM  FORMULATION 


General 


An  infinite  circular  thin-walled  pipe  contains  and  is  also 
surrounded  by  fluids  of  possibly  differing  densities  and  sound  veloci¬ 
ties.  The  displacements  of  the  pipe's  wall  satisfy  an  eighth-order 
shell  theory  t  2  ],  and  the  interior  and  exterior  sound  pressures  satisfy 
the  Helmholtz  wave-equation.  The  pipe  is  excited  by  time-harmonic 
surface  stresses  caused  by  mechanical  point  forces  or  interior  acoustic 
sources.  The  time  factor  exp(-iwt)  is  omitted  from  all  of  the 
equations.  The  geometry  is  shown  in  Figure  1. 
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Pipe  Displacements  and  Fluid  Pressures 


The  axial,  tangential  and  radial  displacements  of  the  pipe’s  wall 
and  the  interior  and  exterior  acoustic  pressures  have  the  following 
Fourier  transform  representations: 
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Formulae  for  the  transforms  U(n,a),  V(n,a),  W(n,a)  and  pa(r,n,a) 
are  to  be  found  in  the  Appendix. 

Far-Field  Acoustic  Pressure 


A  stationary  phase  evaluation  of  the  integral  representation  of 
the  exterior  pressure  yields  the  far-field  pressure  as 


pf (R,8,0) 


W(n,k2cos  B)exp(-imr/2) ,cos(n0) 


-i(DP2c2.exp(ik2R)/irR.  2  . . .  . 

n=0  sin(8)H’ (ak?sin  B) 
n  * 


(2.4) 


Radiated  Power 


The  acoustic  power  radiated  into  the  exterior  fluid  is  defined  as 


P  = 


r* 


•2ir 

[  |p^(R,8,0) |2/p2c2]  R2sin(8)d8d0 


o  Jo 


(2.5) 


where  the  customary  multiplication  factor  of  one-half  has  been  omitted 
because  the  excitation  is  to  be  specified  in  its  rms  form.  Substitu¬ 
ting  equation  (2.4)  into  equation  (2.5)  gives,  after  carrying  out  the 
integration  in  the  0-coordinate, 


P 


(2(i)2p2c2/it)  2 
n-0 


[  |W(n,k2cos  8)  |2/sin(8)  |l-P(ak2sin  8)  |2] 


0 


dS 

(2.6) 


When  the  excitation  is  located  at  za0,  the  symmetry  of  the  sound  field 
about  8air/2  may  be  used  to  halve  the  computational  time. 
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Numerical  Evaluation  of  Integrals 


The  integrals  in  equations  (2. 1-2.3)  and  (2.6)  must  be  evaluated 
numerically  because  closed-form  expressions  are  not  available.  Their 
numerical  approximation,  on  truncating  the  infinite  limits  in  equations 
(2. 1-2. 3)  to  finite  values,  is  based  upon  a  simple  adaptive  quadrature 
scheme  that  has  been  used  successfully  to  evaluate  the  response  and 
acoustic  field  of  a  point-excited  plate  [ 3  ]. 

The  scheme  splits  the  range  of  integration  into  a  selected  number 
of  equal  intervals  each  of  which  being  repeatedly  halved  until  an 
absolute  or  a  relative  convergence  test  is  met  by  successive  approxi¬ 
mations  to  the  integral  in  the  interval  -  Gaussian  quadrature  of  order 
two  is  suitable  for  those  cases  where  low  accuracy  is  sufficient. 

Vogel  &  Feit  1 4  ]  have  also  used  an  adaptive  scheme  to  evaluate  the 
vibration  response  of  a  pipe  with  exterior  fluid  loading  and  point 
force  excitation  -  their  mesh  size  being  controlled  by  the  second 
derivative  of  the  integrand. 

3.  NUMERICAL  EXAMPLES 


General 


Plots  of  wavenumber,  acoustic  power,  acceleration  and  acoustic 
pressure  versus  frequency  are  shown  in  Figures  2-9  for  the  particular 
case  of  a  water-filled  steel  pipe  that  is  surrounded  by  air.  The 
material  and  geometric  constants,  in  SI  units,  used  in  the  computations 
are  as  follows: 

pipe:  E-19.5E10,  cr=0.29,  p=7700.0,  h=2.54E -2,  a=22.2 3E-2,  n=0.02 

water:  p=1000.0,  c=1500.0 

air:  P=1„21,  c=343.0 

Previous  work  ( 1  ]  has  presented  plots  of  far-field  acoustic  radiation 
versus  frequency,  for  the  various  combinations  of  the  fluids,  air  and 
water,  and  the  separate  cases  of  point  force  and  source  excitation. 

Wavenumber  Plots 


Figures  2  and  3  show  the  real  branches  of  the  wavenumber  versus 
frequency  plots  for  the  circumferential  harmonics  n=0-2.  The  loading 
caused  by  the  surrounding  air  was  neglected.  The  frequency  range 
extends  to  1.4  times  the  ’ring'  frequency  of  the  pipe.  The  physical 
interpretation  of  wavenumber-frequency  plots  is  discussed  in  some 
detail  elsewhere  [5]j  here,  only  a  brief  explanation  of  the  'cut-on' 
nature  of  the  waves  is  necessary. 

Figure  2  contains  the  plots  for  a  water-filled  pipe.  First,  for 
n=0,  the  branches  labelled  1-4  are  close  to  a  plane  fluid  wave,  an 
axial  wave  in  the  pipe,  a  radial  flexural  wave  in  the  pipe  and  a  fluid 
wave,  respectively:  the  branch  '0'  is  a  wave  of  pure  torsion. 
Secondly,  for  n=1,  the  branches  are  the  familiar  beam  flexural  wave, 
a  fluid  wave  and  a  pipe  torsional  wave,  respectively.  Finally,  for 
n=2,  the  branches  are  the  radial  shell  flexural  wave,  a  fluid  wave  and 
a  torsional  wave,  respectively. 
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Figure  3  contains  the  wavenumber  versus  frequency  plots  for  the  separate 
cases  of  the  pipe  in-vacuo,  a  rigid-walled  waveguide  of  water  and  a 
pressure  release  waveguide  of  water.  These  plots,  when  compared  with 
Figure  2,  show  the  effect  of  the  interior  water  on  the  real  branches  of 
the  dispersion  relation.  They  are  also  useful  if  the  concept  of  equal 
wavenunber  in  the  fluid  and  pipe  is  of  value  to  the  interpretation  of 
the  responses  of  the  coupled  fluid-pipe  system. 

Airborne  Acoustic  Power 


Figures  4-5  show  the  airborne  acoustic  power  radiated  from  the 
pipe’s  wall  in  the  n=0-2  circumferential  harmonics,  for  the  separate 
cases  of  point-source  and  point-force  excitation.  The  sound  power 
decibel  level  is  defined  as  10.1og(P)  +  120. 

In  Figure  4,  the  excitation  is  a  point  source,  located  at  the 
cylindrical  coordinates  (2a/3,0,0),  whose  free-field  level  in  the 
interior  fluid  is  120  dB  ref  1  micropascal  at  1m.  The  'peaks’  in  the 
spectra  occur  at  the  frequencies  at  which  wave-branches  ’cut-on’.  In 
particular,  a  large  jump  in  the  acoustic  power  occurs  at  the  frequency 
of  the  n=2  harmonic  where  an  essentially  fluid-type  wave  ’cuts-on’. 

At  frequencies  below  the  'cut-on'  frequency  of  the  n=2  flexural  mode, 
the  total  power  is  dominated  by  the  contribution  of  the  n=0  harmonic. 

In  Figure  5,  the  excitation  is  a  radial  1  N  point  force  that  is 
located  at  the  cylindrical  coordinates  (a, 0,0).  Again,  the  peaks  in 
the  spectra  occur  at  frequencies  at  which  the  wave-branches  ’cut-on’. 
However,  the  jump  in  level  at  the  ’cut-on’  frequency  of  the  n=2  fluid- 
type  wave  is  less  severe,  and  the  radiated  power  of  the  n=0  harmonic 
goes  to  zero  at  low  frequencies  where  the  n=1  harmonic  is  dominant. 

The  effect,  on  the  power  radiation,  of  additional  damping  in  the 
pipe’s  wall  was  investigated  by  increasing  the  hysteretic  loss  factor 
from  0.02  to  0.10.  The  power  levels  dropped  by  7  dB  which  indicates 
a  power  reduction  of  10. log ( 112/11) •  It  must  be  borne  in  mind  that  this 
power  reduction  would  not  necessarily  have  been  obtained  for  the  case 
of  a  water-filled  pipe  radiating  into  water,  because  the  radiation  loss 
factor  could  be  of  the  same  order  of  magnitude  as  the  structural  loss 
factor. 

Radial  Acceleration  of  Pipe's  Wall 

The  radial  acceleration  of  the  pipe's  wall  in  the  n=0-2  harmonics 
is  plotted  in  Figures  6-7  for  the  case  of  point  force  excitation  located 
at  (a, 0,0).  The  acceleration  decibel  levels  are  defined  as 

20.1og(|tt|)  +  100. 

Figure  6  gives  the  drive  point  acceleration.  It  is  evident  that 
the  maximum  responses  occur  at  the  'cut-on'  frequencies  of  the  modes. 

The  plots,  on  adjustment  to  velocity  response,  are  of  similar  appearance 
to  the  power  radiation  plots  of  Figure  5. 

Figure  7  gives  the  acceleration  response  at  an  axial  distance  of 
4m  from  the  radial  drive  force.  The  responses  are  now  quite  complex  at 
the  higher  frequencies,  presumably  due  to  modal  interference  effects 
and  the  changing  nature  of  the  wave-branches  with  frequency  [5].  It 
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is  difficult  to  interpret  these  plots  without  the  availability  of  the 
eigenvectors. 

Acoustic  Pressures  on  Pipe’s  Surfaces 

Figures  8-9  show  the  acoustic  pressures  on  the  pipe’s  outer  and 
inner  surfaces,  respectively.  The  excitation  is  a  radial  point  force 
and  the  pressures  are  computed  at  an  axial  distance  of  4m  from  the 
force,  the  decibel  level  being  20,log(Jp|)  +  120. 

The  plots  of  the  pressure  on  the  outside  of  the  pipe.  Figure  8, 
are  of  similar  appearance  to  the  transfer  acceleration.  Figure  7,  At 
the  higher  frequencies,  the  maxima  and  minima  in  the  plots  of  the 
inside  pressure.  Figure  9,  do  not  coincide  in  frequency  with  those  of 
the  outer  pressure  plots.  Indeed,  it  is  curious  that  many  maxima  and 
minima  are  reversed  in  the  two  figures. 
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APPENDIX 


The  Spectral  Displacements 


The  excitation  is  either  a  radial  point  force,  of  magnitude  Fq,  that 

is  located  at  the  cylindrical  coordinate  (a, 0,0),  or  it  is  an  interior 
point  source  of  sound,  of  free-field  pressure  Pgexptik]Rg)/Rg,  that  is 

located  at  (Xg,0,0).  The  'spectral*  displacements  are  obtained,  via 

Novozhilov  shell  theory  [2],  as  the  solution  of  the  matrix  equation 


where 

^11  =  Ej[  a2+n2C1-cr)/2a2  I  -  p  ho2 

s 

s12  =  “Eiian(1+c)/2a 
Si3  =  -Ejiaa/a 
s2l  =  "s12 

522  =  Ei[a2(1-o)/2  +  n2/a2  +  2a2g2(1  -a)  +  (i2n2/a2]  -  pghu)2  (A2) 

523  “  Ej[n/a2  +  na202(2-c)  +  n3£S2/a2  I 

531  =  ~s13 

532  =  S23 

533  *  Eil 1/a2  ♦  a402a2  ♦  2n2o2B2  ♦  n4B2/a2  ]  -  p  w2h 

An  axial  point  force  excitation  is  obtained  by  setting  the  right-hand  side 
of  the  matrix  equation  to  [F0en/2ira,0,0  ]T  . 
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The  Interior  Spectral  Pressure  [  1  ] 

The  interior  ’spectral*  pressure  is  given  by  the  formula 

pjr.n.a)  =  Pjoj2!  Jn(p1  r)/pj  a)  ]  W(n,a)  +  A(n,a)  C  A3) 

where  A(n,a)  is  identically  zero  for  point  force  excitation,  and 

A(n,a)  ■  7T pQen [ Jn (y j x0 ) /J^ Cp j a )  ]  [ 3n (y  1  r) Cp i  a)-Vn (y j rJJ^Cyj a)  ],  r>x^ 

(A4) 

A(n,a)  =  Trp0en[Jn(vi1r)/J^(pia)  ]  [  Jn(u1xQ)V^(yia)-Yn(p1x0)3[!|(y1a)]  ,  r<*g 

CAS) 


for  a  point  source  excitation. 

An  alternative  to  equation  (2.2)  for  the  particular  case  of  point 
source  excitation  is 

-  r 


PiCr.a.z)  =  PgexpCikjRgl/Rg  +  (1/2*)^  cos  (n0) 


b(r,n,a)exp(iaz)  da 


(A6) 


where 

b(r,n,a)  =  [  Jn(y1r)/JUp1a)  ]  [  p^o^WCn.aJ/yj  -  uip^e^Cy^^fT  (y^)) 

(A7) 
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(J  )  x  =  R  .  si  n  (  0)  cost  0  ) 

(f  )  y  a  Rsin(0)sin(jJ) 

z  =  R . cos(0  ) 


D  COORDINATE  SYSTEMS 
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PRESSURE  ON  OUTSIDE  OF  PIPE.  Z=4.8H 
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